Domain excitations in spin-Peierls systems. 
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We study a model of a Spin-Peierls material consisting 
of a set of antiferromagnetic Heisenberg chains coupled with 
phonons and interacting among them via an inter-chain elastic 
coupling. The excitation spectrum is analyzed by bosoniza- 
tion techniques and the self-harmonic approximation. The 
elementary excitation is the creation of a localized domain 
structure where the dimerized order is the opposite to the 
one of the surroundings. It is a triplet excitation whose for- 
mation energy is smaller than the magnon gap. Magnetic 
internal excitations of the domain are possible and give the 
further excitations of the system. We discuss these results 
in the context of recent experimental measurements on the 
inorganic Spin-Peierls compound CuGeOa 



The recent discovery of the first inorganic Spin- 
Peierls(SP) compound CuGeOs has renewed the interest 
in the subject, of quasi-one-dimensional spin-phonon cou- 
pled systemsEl. In spite of the intense experimental and 
theoretical activity devoted to the study of this system 
neither the excitations in the low temperature phase nor 
the mechanism of the Spin-Peierls transition itself are 
well understood yet. Most studies have used as a model 
Hamiltonian a Heisenberg antiferromagnetic chain with 
alternating^aoupling and next nearest neighbors (nnn) 
InteractionBu. The Spin-Peierls transition is supposed to 
arise from a competition between the elastic energy cost 
necessary to dimerize the chain and the magnetic energy 
gain in this process. The dynamics of the phonons are 
then supposed to be independent of the magnetic subsys- 
tem and an extreme adiabatic approximation is assumed. 

In spite of the success to reproduce both thermody- 
namical as well as dynamical properties of Spin-Peierls 
systems by this approach, there are now several indica- 
tion that the Spin-lattice interaction in CuGeOs is strong 
enough to make insufficient this approximation. No sc 
phononic mode related to dimerization has been founc 
This fact point toward an order-disorder type transition 
where nonlinear excitations (domain wall) different from 
phonons drive the structural transition. Otherwise, the 
phase diagram of the system in presence of a magnetic 
field is now well established. Above a critical value of the 
magnetic field the system undergoes a transition from 



the uniform dimerized into an incommensurate phase. 
The incommensurate lattice pattern has been measured 
by X-ray_experiments and interpreted as soliton lattice 
structureB. The role of solitons in CuGeOa was indeed 
previously emphasizedtl. They could give an unified pic- 
ture of the above features. They could also be relevant 
to explain the rapid reduction of the SP temperature 
by dooipg and the apparition of an antiferromagnetic 
phase&Q. 

It seems then natural to build the excitation spectrum 
of SP systems on the basis of some kind of solitonic exci- 
tation (in this work we use the term soliton in the general 
sense of a finite energy localized configuration which is 
simultaneously magnetic and structural excitations, we 
call kinks the ID solitons). The present work gives some 
insight in this direction. 

Moreover, recent neutron scattering studies have given 
a detailed information about magnetic excitations in 
CuGeOaO. They show a 'double gap' structure, in ad- 
dition to a dispersive triplet excitation there is another 
gap which separates the first peaks from the band edge 
of a continuum. This continuum has been interpreted 
as due to kink excitationsEl. It was stated that ao free 
kinks could exist due to the inter-chain couplingtl. The 
inter-chain coupling were taken into account in previous 
works as nproviding a linear confined potential between 
the kinksH. This potential arises from a mean field ap- 
proximation to the inter-chain coupling. Particularly in 
a recent work a laddeii.of bound kink-antikink states was 
shown for this modellld. 

In this paper we study the formation of solitonic struc- 
tures in SP system including the three dimensional char- 
acter of the phonon field. We use bosonization techniques 
to account for the low energy magnetic excitations. We 
show that the excitation gap of the system corresponds 
to create a domain-like localized structure in which the 
dimerization is in antiphase to the one in the bulk mate- 
rial. No independent kinks are possible. Instead of that 
the walls of the domain always find an equilibrium situa- 
tion depending on the value of the inter-chain interaction. 
Internal excited states of the domain are the higher ex- 
cited levels of the system. They could be connected with 
the continuum seen in CuGeOa. 

Let us consider a system of Heisenberg antiferromag- 
netic chains immersed in the phonon field of the mate- 



rials. The spin-phonon coupling arises from the modu- 
lation of the magnetic exchange by the lattice motion. 
Therefore we will focus on the following spin-phonon 
Hamiltonian: 



H — Hph + H, 
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we denote by i the site in a j-th chain. S^ are spin-1/2 
operators of the i,j ion. Hph represents our simplified 
model for the phonons. It contains a scalar coordinate 
ul and its conjugate momentum P^ which are supposed 
to be the relevant ion coordinates for the dimerization 
process. Notice that, in CuGeOa, these are mainly re- 
lated with the displacement of the O ions. K\^ and K± 
are the longitudinal and transversal elastic couplings. We 
have included only the coupling in one direction, gener- 
alization to the other direction is straightforward and 
produces only a redefinition of the quantity B(j) of Eq. 
(O). This transversal coupling gives the three dimen- 
sional character to the phonons field which has been jaa- 
glected in previous studies of solitons in SP systemllil. 
It is essential for our considerations. The model Hamil- 
tonian given by (m) does not include magnetic interac- 
tion between the chains. We will discuss latter on the 
possible effect of this interaction on the excitation spec- 
trum. Moreover we have not introduced a nnn interac- 
tion between the spins of a chain. Even though, it play 
an important role in reproducing several experimental 
properties in CuGeOa we do not expect any qualitative 
change of the results presented below which are based on 
the solitonic excitatiops-smd not strongly affected by this 
frustrated interactiorOliil. 

A similar model of the one given by Eq. [^) were re- 



cently studied by Essler, Tsvelik and DelfinoLJ using rig- 



orous field theoretical results 
of ref . 
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However, as the authors 
have integrated out the lattice coordinates and 
neglected the dynamic of the phonons they did not take 
into account the mixed excitation we will discuss in the 
present paper. We remark that in CuGeOa the scale of 
magnetic and phononic excitation are of the same or- 
der of magnitude suggesting that the elimination of the 
phonons degree of freedom is not suitable. 

Owing to the one-dimensional character for the mag- 
netic interaction in (|l|) we can resort to the bosonization 
techniques to tate into account the low energy excitations 
of the systemalj. The spin variables are approximately 
represented by the bosonic fields (f>^ (x) and its conjugated 
momentum H-' (x). In addition we retain only the phonon 



modes producing a smooth deviation of the dimerized 
pattern. So, we make the replacement (— 1)*m^ —^ u^{x). 
The effective Hamiltonian {H^ff = Hph + Hmg) is given 
by: 

2 
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Vs is the spin wave iielocity and rj the exponent of the 
correlation functiongl3'tll. For the isotropic Heisenberg 
model with nn interaction we have Vs = 7r/2 and rj — 2. 
We analyze the excitations of -ffe// by using a 
self-consistent harmonic approximation(SHA) whiohihas 
been shown to be reliable for Spin-Peierls systemstll. It 
is nearly connected to the semi-classical i»ethods widely 
used in the study of nonlinear field theoryt3. In the SHA 
the fields cj)^ are split into a classical variable ^q ^^'^ its 
quantum fluctuation (j)^ . Terms up to quadratic order in 
(j)i are retained and the annulment of the first order term 
is required . The resulting equations are: 



— {dl<l>'o) + w^e-<*">/2 ^^g^, ^ 
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(0J ) is the ground state expectation value and the last 
equation is the classical equation for the lattice coordi- 
nates u^ . Uniform solutions of (g) and (|^) correspond to 
the homogeneous dimerized configuration of the chains. 
They are given by: 



-{mod.2T:) , u-' (x) = uq 
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the transverse coupling is not effective in this case and 
the state is the same as the one of a single chain problem. 
Excitations over this state are given by: 



u;{k) = Us(^ +fco) 



2N1/2 



(9) 



they represent a band of triplet excitations separated 
from the ground state by the gap A — v^kQ — (ttuq)^^'^ . 
These magnon-like excitations are usually characterized 
as the elementary excitations of a Spin-Peierls system. 

However, there is a lower excitation energy of the sys- 
tem which arises when inhomogeneous solutions of (H) 
and (|^) are allowed. To take it into account we invert 



eq. (0) to obtain u^x) and then replace it in eq (g). The 
result is: 



u^ (x) = "0 ^ B(/ - j) sin 0i?, (x) 
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is the relative inter-chain elastic couplings and 

£, — Vs/A will become the characteristic width of the 
domain wall(see below). 

Eq. (11^) and (lO) are our basic system of differential 
equations to search solitonic excitations of the system 
(each equation is labeled by the index j of the chain). 
They should be solutions of the system under the re- 
quirement of total finite energy respect to the uniform 
configuration. The inter-chain elastic coupling strongly 
limits the form of these solutions. For example, free kinks 
do not fulfill this condition because they create an infi- 
nite zone where the transverse couplings are active. The 
fields u^{x) and cj)^ {x) should go as a; ^- ±oo to the values 
they have in the uniform configuration (0) . The simplest 
solution satisfying these conditions can be constructed as 
follows. We freeze the value of (J>q in all the chains but 
the 0-th to its value in the homogeneous SP state, i.e. 
(J)q = n/2 for j 7^ 0. The system ( pT| ) then reduce to a 
simple equation over the 0-chain. It corresponds to the 
classical equation for an effective one-chain problem with 
renormalized parameters and an additional term favoring 
the dimerization phase of the neighbors chains. We solve 
the equation by direct integration imposing the previ- 
ously discussed boundary conditions. The result is: 

t{x) = sin((/)g) = 1 — 2 cosh^[a;o/C] sech[(a:: — xq)/^] 



could not collapse this value never reaches the double of 
the width of the wall(2^). Note that for each inter-chain 
coupling we have an equilibrium distance. We do not 
find a linear confinement poteiitia|l between the kinks as 
it is assumed in other analysia§ll3 of this problem. The 
domains turn out to be the elementary excitations of the 
system. They will move as a whole when the translational 
invariance be restore. This is a triplet excitation because 
its total Sz could be or ±1. Note that in the bosonic 
representation the total magnetization is given by 5*°* — 

2^V.(,/,^(00)-<^-(-^)). 
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FIG. 1. (a)function t{x) defined in (|13[). It gives the mag- 
netic profile of the domain excitation. The insertion shows 
the e-dependence of the parameter 2xq (related to the do- 
main's size) and & named the distance between the zeros of 
t(x). b)Lattice pattern corresponds to ( [f5| ) for the e — 0.57. 
We show the deformation of the 0-chain and its neighbours 
chains. 

The lattice deformation pattern is obtained from (ll 



v?{x) = uoB(j)sin0SJ + (1 - B(j))7/o 



(15) 



sech[(2; + a;o)/C] (13) 
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Xo 



log 



1 + B(0) + 2v/B(0) 
l-B(O) 



B(o) = xj^/YTl 
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This solution is shown in fig. (|l|). Wc call it a domain. 
2xo measures the size of the domain configuration ( see 
the insertion of fig. (|l|)). It goes to infinity as the chain 
get decoupled (e — > 0). In this limit, expression ( |l^ ) 
becomes a product of two tanh-kink form. A small e 
produces a rapid accommodation of the walls one near 
the other in order to reduce the inter-chain energy. For 
a greater e, 2xo departs from the distance between the 
zeros of t{x) called 8 in the inserted figure. As the walls 



We show this configuration in fig. (|i|b). The displace- 
ment amplitude in the intermediate zone of the 0-th chain 
is smaller than uq- The distortion is not restricted to the 
0-chain as for the bosonic field 0o • The neighbors chains 
are slightly distorted respect to the uniform dimerized 
pattern. This distortion is indeed very weak as-the re- 
sult of the rapid fall of the value of B(j) with jll3. This 
behavior together with the asymmetry seen between the 
magnetic and lattice pattern is possible as a consequence 
of neglecting the inter-chain magnetic interaction. We 
should expect that inclusion of this interaction will smear 
out the deformation producing a kind of two (or three) 
dimensional domain with dimerization pattern opposite 
to the one of the surrounding material. Note that in our 
model spin-one excitations live in a chain but on turning 
off the transversal magnetic exchange a triplet localized 
state could appear exciting simultaneously several chains. 
Summarizing, a general solution of (O) in presence of a 



magnetic interaction will be an enlarged domain in the 
direction of the chains. 

We now turn to the study of the quantum states of 
the theory in terms of our classical solution. The first 
contribution comes from the creation energy of the do- 
main. It is the difference of the values of the classical 
energy plus the zero point energy of the fields obtained 
with and without the domain. The classical energy is 
evaluated by replacing ( p^ ) and ( p^ ) in (|^) and ^ and 
subtracting the energy of the uniform configuration. Af- 
ter a long but straightforward calculation we obtain: 
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The e-dependence of this classical energy in shown in fig. 
(ga). For e ^ we have Ed = 2.Jk '^hich is twice 
the classical creation energy of a kink in the single chain 
problem (expressions (4.10) and (4.11) of ref. 11 with the 
appropriate redefinitions of the constants). When the 
inter-chain coupling is switch on, the walls get closer and 
its total energy shghtly increases. Note that in fig.(||a) we 
are showing the energies in units of the magnon gap A, 
therefore for any reasonable choice of e it will be more 
favorable to create a domain rather than to a excite a 
magnon. This result is at the hearth of our claim that the 
domains are the lower energy excitations of SP systems 
and will survive when quantum fluctuation are included. 




(b) 



FIG. 2. (a) Differents contribution to the domain energy 
(Ed)- Ed classical energy, Eb3 contribution of the bound state 
of the potential (|l7|), Econt contribution of the continuus(see 
text), (b) Energy levels diagram for the parameter fitted for 
CuGeOa. Dashed lines denote the excited levels of the domain 



We now go to the inclusion of these quantum fluctu- 
ations. The fluctuation operator is the quadratic form 
in the fields (fP of the SHA Hamiltonian. Its eigenvalue 
problem corresponds to a set of Schroedinger-like equa- 
tions for the following potentials: 



V°{x) = {'B{QMx)-l) + l)t{x) 



V\x)=^{j){t{x)-l) + l 



(17) 



(18) 



Once the eigenvalues {Ex) are known the frequency oscil- 
lations in presence of the domain could be computed as 
'^d(A) =_Ay/"E\. The quantum correction to the classical 



energy (16) is the difference between the sum of the zero 
point energies of these oscillators and the ones given by 
(0) i.e. in absence of the domain (the cutoff dependence 

is eliminated when the tppi containing {(j)i ) in the SHA 
Hamiltonian is includedllij) . 

Equation (O) defines a quantum mechanics problem 
of a particle in a double well potential, whereas, Eq. 
(O) the one of weak wells. We have solved numerically 
the eigenvalues problem. The bound states has been ob- 
tained by means of a nodes counting algorithm. Numer- 
ical details of the method will be given elsewhere. For 
the double well potential (p^), we have found two bound 
states (called Ei and E2)- For e — > 0, the wells get far 
apart and the two levels collapse in a degenerate bound 
state with the same-energy as the bound state of the 
single kink problemO. The larger e, the bigger is the 
splitting due to the tunneling between the wells. 

Potentials (^8|) have only one bound state near the bor- 
der of the continuum due to the weakness of the potential. 
We neglect its contribution to the total creation energy 
as well as the slightly distortion of the continuum. We 
will discuss later on the relevance of these bound states 
as excited levels of the domain. 

The contribution of the continuum levels of ( [l7|) could 
be computed once the phase shift 6{k) of these states 
is known. We calculate its by numerical integration of 
the equation starting from the origin and matching its 
behavior to its asymptotic behavior at x — > 00. The final 
expression for the total contribution of these states is: 
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(19) 



The last integral has been numerically evaluated. In fig. 
(^) we show the e-evolution of the different pieces con- 
tributing to the total creation energy (iJ^). The main de- 
pendence comes from the classical energy. We find, again, 
that in the limit of decoupled chains our excitation en- 
ergy corresponds to the total energy of two independent 
kinksB. 

What about CuGeOs? As we have stated in the intro- 
duction, neutron scattering experiments show a clear dis- 
persive excitation at 2.1 meV. This excitation was previ- 
ously analyzed on the basis of a dimcrized and frustrated 



Heisenberg chainB'Olia. Taking into account the previ- 
ous results we associate the energy gap with the energy 
formation of a domain. To fix parameters we use the 
experimentally fitted values J = 120K and uo — 0.042. 
With these values we have A = 2.68 meV larger than the 
gap. We can therefore choose e — 0.57 giving Ed equal 
to the measure gap. For this e the walls are near and 
the domain seems as a local depression of the SP order 
(see fig.0). Inclusion of nnn interaction will change the 
parameters. The next energy levels of the system are 
internal excitations of the domain. They correspond to 
the bound states previously discussed. Fig. (||b) shows a 
level diagram for the low energies excitation of the sys- 
tem. It is natural to associate the transitions to the ex- 
cited levels with the continuum seen in CuGeOs. We 
recall that switching on the inter-chain magnetic cou- 
pling will smoother of the domain deformation in the 
direction perpendicular to the chains. Therefore we ex- 
pect that the internal excitations will form a continuous 
band. Note also that a two domain continuum is possi- 
ble starting near 2E4. To give precise predictions of the 
effect of the domain formation on the neutron scattering 
spectra the detailed information about the dynamics cor- 
relation functions will be needed. Note in addition that 
as our domain is referred to a given state of dimerization 
degeneration between the states at g = and g = tt is 
expected. 

Let us assume a level excitations diagram for CuGeOa 
as the one of fig. (Ep) . As a finite density of domain will 
be thermally created our model predict transitions cor- 
responding to rise the domain to its excited levels. The 
first one should appear as a low energy peak at about 
u)d{^) — 1-78 meV in the optical response of the system 
and to increase its intensity with the temperatur|e-|Up to 
Tsp- Recent inelastic light scattering experienceslij have 
identify a new resonance in agreement with this predic- 
tion. However, note that this transition were assigned 
in Ref 03 to a three-magnon process on the externally 
dimerized chain. More experimental and theoretical work 
will be need to decide between these interpretations. 

Finally, we have discussed in this paper only the mag- 
netic excitations of the system. However the domain is 
a mixed state between the spins and phonons. So we 
expect that the phonon response will be sensitive to the 
domain formations. The apparition of shoulder or satel- 
lite peaks related to the frequency oscillations of the ions 
around the domain should be expectedc2l. 

In summary, we have identified a new excitation mode 
of Spin-Peierls systems. It gives a natural extension of 
the usual singlet-triplet excitation of the statically dimer- 
ized chain when the lattice relaxation is allowed. As the 
excited spins are not confined to a dimer, internal excita- 
tion are possible. We claim that these excitations could 
account for some of the features seen in recent experi- 
ments. 
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